ASYMPTOTIC ANALYSIS IN A GAS-SOLID COMBUSTION 
MODEL WITH PATTERN FORMATION 



CLAUDE-MICHEL BRAUNER, LINA HU, AND L. LORENZI 

Abstract. We consider a free interface problem which stems from a solid-gas 
model in combustion with pattern formation. We derive a third-order, fully 
nonlinear, self-consistent equation for the flame front. Asymptotic methods 
reveal that the interface approaches a solution of the Kuramoto-Sivashinsky 
equation. Numerical results are presented which illustrate the dynamics. 



1. Introduction 

Combustion phenomena are particularly important for Science and Industry, as 
J.-L. Lions pointed it out in his foreword to the special issue of the CNRS "Images 
des Mathematiques" in 1996 (PQ). Flames constitute a complex physical system 
involving fluid dynamics and multistep chemical kinetics (see e.g., [9])- In the mid- 
dle of the 20th century, the Russian School, which included Frank-Kamenetskii 
and Zel'dovich, used formal asymptotics based on large activation energy to write 
simpler descriptions of such a reactive system. Later, the development of system- 
atic asymptotic techniques during the 1960s opened the way towards revealing an 
underlying simplicity in many combustion processes. Eventually, the full power 
of asymptotical analysis has been realized by modern singular perturbation the- 
ory. J.-L. Lions was the first one to put these formalities on a rigorous basis in his 
seminal monograph "Perturbations singulieres dans les problemes aux limites et en 
controle optimal" (|14j). 

In short, the small perturbation parameter in activation-energy asymptotics is 
the inverse of the normalized activation energy, the Zel'dovich number (3, In the 
limit j3 +oo, the flame front reduces to a free interface. The laminar flames of 
low-Lewis-number premixtures are known to display diffusive-thermal instability 
responsible for the formation of a non-steady cellular structure (see |23]). when 
the Lewis number Le (the ratio of thermal and molecular diffusivities) is such that 
Le < 1. From an asymptotical viewpoint, one combines the limit of large activation 
energy with the requirement that a = |/3(1 — Le) remains bounded: in the Near 
Equidiffusive Flame model (NEF), and 1 — Le are asymptotically of the same 
order of magnitude (see [2T]V 

A very challenging problem is the derivation of a single equation for the (per- 
turbation of the) free interface, which may capture most of the dynamics and, 
as a consequence, yields a reduction of the effective dimensionality of the system. 
Asymptotical methods are again the main tool: in a set of conveniently rescaled 
dependent and independent variables, the flame front is asymptotically represented 
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(see [22.) by a solution of the Kuramoto-Sivashinsky equation: 



$ r + 4$ 



wivv 



(K-S) 



This equation has received considerable attention from the mathematical commu- 
nity (see [H]), especially for its ability to generate a cellular structure, pattern 
formation, and chaotic behavior in appropriate range of parameters (see |13|). We 
refer to 0-[7] for a rigorous mathematical approach to the derivation of (K-S). 

In this paper, we consider a model in gas-solid combustion, proposed in |12j . 
This model was motivated by experimental studies of Zik and Moses (see [25]). 
who observed a striking fingering pattern in flame spread over thin solid fuels. The 
phenomenon was interpreted in terms of the diffusive instability similar to that 
occurring in laminar flames of low-Lewis-number premixtures. As we show below, 
the gas-solid and premixed gas systems share some common asymptotic features, 
especially K-S equation. 

The free interface system for the scaled temperature 9, the excess enthalpy S, 
the prescribed flow intensity U (with < U < 1), and the moving front x = £(i, y), 
reads: 



u» - A,, 

ox 
9 = 1, 

89 u dS_ 
dt dx 



x > (,(t,y), 
= AS-aA9, 



(1.1) 
(1.2) 

(1.3) 



System Ql.l|) - (|1.3|l is coupled with the following jump conditions for the normal 
derivatives of 9 and S: 



89 



8n 



■exp(S), 



dS 



dn 



a 



89 
dn 



(1.4) 



It is not difficult to show that System (jl.ip ~ p.4p admits a planar Traveling 
Wave solution, with velocity —V = U\nU. Setting x' = x + Vt, the Traveling 
Wave solution reads: 

{ exp(Ux'), x' < 0, 
9(x') = I ' 

I 1, x'>0, 



S(x') 



(a - In U)Ux' exp(Ux') + (In U) exp(Ux'), x' < 0, 



[ bxU, x' > 0. 

As usual one fixes the moving front: we set 

S(t,y) = -Vt + <p(t,y), x' = x-Z(t,y), 

where ip is the perturbation of the front. In this new framework the system (jl.ll) - 
(fL"5j) reads: 



U9 X , = A v 9, x 1 < 0, 
9 = 1, x' > 0, 

9 t + (V- y t )9 x , + US X , = A V S - aA v 9, x' f 0, 



(1.5) 
(1.6) 
(1.7) 



where 



A,p = (1 + {<-Py) 2 )D X > X l + Dyy — iPyyD X l - 2if y D X > y . 
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The front is now fixed at x' — 0. The first jump condition in (| 1 .4[) reads: 

89' 



1 



the second one: 



(<Pv) 2 



dS 



dx' 



dx' 



= -exp(S), 



de 



dx' 



(1.8) 



(1.9) 



We will consider a quasi-steady version of the model, motivated by the fact 
that, in similar problems, not far from the instability threshold, the respective time 
derivatives (if any) of the temperature and enthalpy exhibit a relatively small effect 
on the solution: the dynamics appears to be essentially driven by the front. We 
can thus introduce a quasi-steady model replacing (|1.5|l - (|1.7|l by 



U9 X > = A p 0, x < 0, 
9 = 1, x'>0, 

(V- Vt )9 x ,+US x ,=A v S-aA v 9, x'^0. 

Next we consider the perturbations of temperature u and enthalpy v. 

9 = 9 + u, S=S + v, 

and, for simplicity, in the equations satisfied by u and v and tp, we keep only linear 
and second-order terms for tp, and first-order terms for u and v. Writing x instead 
of x' , to avoid cumbersome notation, some (easy) computations reveal that the 
triplet (u, v,tp) solves the differential equations 

Uu x - Au = {A v - A)9, x < 0, 

Vu x - A{v ~ au) + Uv x - ip t 8 x = (A v - A)(S - ad), x ^ 0, 
where u = in [0, +00) and 

(Ap - A)9 = {U(ip y ) 2 - <p yy )Ue Ux , 

(A<p - A)(S - aff) 

( {tp y ) 2 {a - In U)U 2 {\ + Ux)e Ux - ip yy (a - In U)U 2 xe Ux , 



0. 



x < 0, 
x > 0. 



The previous system is endowed with a set of boundary conditions. First, the 
continuity of 9 and S at the front yields the equation 



u(0~) 



0. 



(recall that u(x) = for x > 0). Second, condition (|1.8[) gives, up to the second- 
order, 



-U+[u x ] = -{l + (^ y ) 2 )-iUe 



v(0) 



1 



- )U l + «(0) + -(«(0)) 



and keeping only the first-order for v we get the condition 

-u x {0-) + Uv(0) = ^{tfyfU. 

Finally, the condition [S x ] — a[9 x ] yields 

[v x ] = -au x {0~). 
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Summing up, the final system is the following one: 

Uu x - Ait = (A v - A)0, x < 0, 

Vu x - A(v - au) + Uv x - cp t 9 x = (A v - A)(5 - d§), x ^ 0, 

u(o~) = H = o, (i.io) 

tf«(0) -u x (Q-) = \{<p v ?U, 
[v x ] = -otu x (0~). 

Throughout this paper we will also use the very convenient notation: 

7 = a — In U. 

First, our goal is to derive a self-consistent equation for the front <p: 

<P t = £/(<p)+^((<fy) 2 ), (1.11) 

where is a third-order, pseudo-differential operator, in contrast to the NEF 
model in gaseous combustion, where the corresponding linear operator is of the 
second-order (see [5]). Another important feature is that the nonlinear term is also 
of the third-order, which means that the equation (jl.lip is fully nonlinear. Here the 
spatial domain is a two-dimensional strip M x [—1/2,1/2], with periodic boundary 
conditions at ±£/2. 

Second, we define a small parameter e = 7 — 1. The main result of this paper 
states in which precise sense the front ip approaches a solution of the Kuramoto- 
Sivashinsky equation when e — > 0: 

Main Theorem Let $0 € H m (— &q/2, £q/2) be a periodic function of period £0. 
Further, let $ be the periodic solution of (|K-S[) (with period £q) on a fixed time 
interval [0,T], satisfying the initial condition ^(O, ■) = $o- If m is large enough, 
then there exists eo = £q(T) G (0,1) such that, for < e < Eq, Equation (|l.lip 
admits a unique classical solution ip on [0, e ^ u2 ] , which is periodic with period 
with respect to y, and satisfies 

v(0,y) = eU- 1 <P o (yV^U), \y\ < '" 



2^feW 

Moreover, there exists a positive constant C such that 
\<p(t,y)-eU- 1 1>(te 2 U 2 ,y^U)\<Cs 2 1 0<t<^-,, \y\ < ' >: 



e 2 U 27 m " 2y/eU' 
for any e € (0,e ]. 

The paper is organized as follows: in Section[2j we proceed to a formal Ansatz in 
the spirit of |22) . defining the rescaled variable ip = Uip/s and expanding ip — tp + 
eip 1 + . . .. It transpires that ip° verifies (IK-S|) . thanks to an elementary solvability 
condition. 

Section [3] is devoted to the derivation of (II. lip , via an explicit computation in 
discrete Fourier variable. The asymptotic analysis in the rescaled variables t = 
T/£ 2 U 2 ,y = r)/y/eU is performed in Section|4l Since the perturbation in (jl.lll) is 
singular as e — > 0, we turn to the equivalent (at fixed e > 0) fourth-order, fully 
nonlinear equation (|1.12p . whose prima facie limit as e — > is Equation (IK-SI) : 

(y/l - ip = -4D vvvv ip - D vv ip 



Ot 



ip) 2 . 



J {(/ - ±eD m )i - 3(1 - 4eD m ) - 4(1 + e) (yj - AeD m - ij } (D v ip) 

(1.12) 
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We prove a priori estimates which constitute the key tool to prove the Main The- 
orem. Finally, numerical computations which illustrate the dynamics in Equation 
(|1.12p are presented in Section [5] 

The stability issue in System (jl.ip - (|1.4[) will be addressed in a forthcoming paper, 
using the methods of [5] and [15]- [15]. 

Notation. Given a (smooth enough) function / : (—1/2,1/2) — > C, we denote by 
/(fc) its fc-th Fourier coefficient, that is, we write 



/(y) = E/(*> fe ^)' ve(-t/2,£/2), 



fc=0 

where {wk} is a complete set of (complex valued) eigenfunctions of the operator 
D yy : H 2 (-l/2,l/2) -> L 2 (-£/2,£/2), 

whose eigenvalues 0,-^,-^,-i|f^,-^,-^,... we label as = -A (^) > 
—Xi(£) = — A 2 (^) > — A 3 (i?) = — \4(£) > . . . Typically, when no confusion may arise, 
we simply write A& instead of Xk (I) . 

For any s > we denote by i?| the usual Sobolev space of order s consisting of 
^-periodic (generalized) functions, i.e., 

{+oo +oo ~j 

fc=0 fc=0 J 

For s — 0, we simply write L 2 instead of and we denote by | • | 2 the usual 
L -norm. 

By the notation /(x, fc) we mean the fc-th Fourier coefficient of the function 
f(x, •). A similar notation is used for functions which depend also on the time 
variable. 

2. A FORMAL ANSATZ 

The aim of this section is to use a formal asymptotic expansion method, in the 
spirit of |22j . The small perturbation parameter e > is defined by: 

a= l + lnf/ + e, i.e., 7 = l + e. (2.1) 

Accordingly, we now introduce scaled dependent and independent variables: 

*=^a> y= ^U' ^ = ^> u = e 2 u u v = e 2 Vl , (2.2) 

and the Ansatz: 

Mi = it? + eu\ + . . . , vi=v°+ev{ + ..., ip = ip° + sip 1 + . . . . 

It is easy to rewrite System (jl.lOp in terms of the rescaled variables. At the zeroth 
order, it comes: 

{ U(u\) x - (u\) xx = -U 2 e u ^° vv , x < 0, 

V{u\) x - {vf) xx + {u\) xx + (In U)(u<i) xx + U(v^) x = -U 3 xe Ux ^ n , x < 0, 

u? = 0, x > 0, 

(«?)» - - 0, x > 0. 

(2.3) 

At x = 0, the following conditions should be satisfied: 

u°(0) = = 0, (2.4a) 
(uS) s (0) - Uvl (0) = 0, (2.4b) 
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[(«?)*] =-(l+lni7)K) x (0). (2.4c) 

We assume that the functions x e~ Ux l 2 u\(x) and x H> e - ' 71 / 2 -^^) are bounded 
in (—00, 0) and in M, respectively. Note that System (12.31) coupled with conditions 
(I2.4a|) and (|2.4bl) is uniquely solvable in the unknowns (ui,Vi), taking ip° as a 
parameter. It turns out that 

u\ = Uxe Ux i>° m , x < 0, 

w° = e Ux iP° v + U(ln U)xe Ux ^ + U 2 x 2 e Ux ^ x < 0, 
u? =0, x > 0, 

*>0. 

One might be tempted to use condition (|2.4c[) to determine function Unfortu- 
nately, whichever ip° is, the triplet (v,®, v®, satisfies this condition. As a matter 
of fact, we are not able to determine uniquely a solution to Problem (|2.3[) - (l2.4cp . 
This situation is not surprising at all in singular perturbation theory, see [10] , |14] . 
To determine ip° one needs to consider the (linear) problem for the first-order terms 
in the asymptotic expansion of u\, v\ and ip. As we show in a while, this problem 
provides a solvability condition, which is just the missing equation for ip° . 
The system for (u\,v\,ip l ) is the following one: 

U(u\) x - (u\) xx - U 2 {u\) m = (C/(^°) 2 - U^ v )Ue Ux , x < 0, 

V{u\) x - {v\) xx - C/ 2 (w?) w + («?) aa + {l + \nU){{u\) xx + t/ 2 K) w ) + U{v\) x 
= U 2 ^ T e Ux + (^) 2 U 2 e Ux + (^) 2 U 3 xe Ux - U 3 ^ v xe Ux - U 3 ^ v xe Ux , 

x < 0, 



{v\) xx + U 2 {vl) m - U{v\) x = 0, x > 0, 

u\ =0, x > 0, 

«i(0)= [vl] =0, 

[(«!)«] = -(l + ln£7)(«l) s (Q)-(ti2)*(0). 

(2.5) 

As above, we assume that the functions x H> e Ux ' 2 u\(x) and x i— > e t7 ' a! ' 2 t>J (a;) are 
bounded in (— oo, 0) and in K, respectively. Using these conditions one can easily 
show that the more general solutions (u\,v\,ip l ) to the differential equations and 
the first boundary condition in (|2.5p are given by 



Uxe Ux ^° vvvv - (^°) 2 + - \u 2 x 2 e Ux ^ mm , x < 0, 



where 



v{ = vl(0)e Ux + Axe Ux + Bx 2 e Ux + Cx 3 e Ux , x<0, 
u\ = 0, x > 0, 
vl=vl(0) + Ux^ mri , x>0, 

A =U(lnU)(^ n - (V^) 2 + - UiP° T - U(^) 2 - 3Ety° 



>l>l>l>l ' 



B =U 2 ^ m + C/ 2 ^ r) - U 2 (^) 2 - \u\\nU)^ mm + \u 2 ^ mv 
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c = - -^u 3 ^ vriv , 

and v\(0) is an arbitrary parameter. Hence, depends on ip . To determine 

both-0 1 and v\(0) we use the last two boundary conditions, which give, respectively, 

-EtyJU, + U(^) 2 - Utfn + Uv{(0) = |«) 2 (2.6) 

and 

Utf n - Uv\ (0) = - Utf - U^ m - W^ mm , (2.7) 

Obviously (|2.6[) and (|2 . T|) is a linear system for (u*(0), ^L) with solvability condi- 
tion: 

^ + V4 + 4^ w + i(v^) 2 = o, 

i.e., '0° verifies a K-S equation. Hence, the Kuramoto-Sivashinsky equation is the 
missing equation, at the zeroth-order, needed to uniquely determine {u\, v®, 

3. A THIRD-ORDER FULLY NONLINEAR PSEUDO-DIFFERENTIAL EQUATION FOR 

THE FRONT 

The aim of this section is the derivation of a self-consistent pseudo-differential 
equation for the front tp. We consider System f|l . 10|) . namely: 

Uu x - Au = Ue x (U(ip y ) 2 - (p yy ), x<0 : 

Vu x - A(v - au) + Uv x = Uip t e x + (a - In U)U 2 (1 + Ux)e Ux (ip y ) 2 

-U 2 {a - \nU)xe Ux (p yy , x < 0, 

Uv x -Av = 0, x>0, 

t*(0) = [v] = 0, 

Uv(0)-u x (0) = \U{ip v ) 2 , 

. [vx] = -au x (0), 

(3.1) 

in a two-dimensional strip K x [—£/2,£/2], with periodicity in the y variable. 

3.1. Computations in discrete Fourier variable. Throughout this section, 
(u, v, ip) is a sufficiently smooth solution of System (I3.1[) such that the functions 

(x, y) h> e~ Ux/2 u{t, x, y), {x, y) H- e~ Ux ^v{t, x, y) 

are bounded in (-oo,0] x [-1/2,1/2] and infix [-£/2,£/2], respectively. 
We start from the first equation in (|3.1|) . namely: 

Uu x -Au= (U(tp y ) 2 -tp yy )Ue x , (3.2) 

and the boundary condition «(•, 0, •) = 0. Applying the Fourier transform to both 
the sides of (|3.2|) we end up with the infinitely many equations 

Uu x (t,x,k) -u xx {t,x,k) + \ k u{t,x,k) = (u '(tp y ) 2 (t,k) + \ k 0(t,kfj Ue Ux , 

for k > 0. For notational convenience we set z/ fe = ~ + \\/U 2 + 4Afe for any k > 0. 

Since u vanishes at x = and tends to as x — > — oo not slower than e Ux ^ 2 , the 
modes «(•,•> k) should enjoy the same properties. Easy computations reveal that 

u(t,x,0) = -UlJpyf{t,0)xe Ux , x<0, 



u(t, x, k) = UiXk)- 1 (U(ip y ) 2 (t, k) + X k 0(t, fc)) (e Ux - e v " x ) 



x < 0, k > 1. 
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Applying the same arguments to the equation for v jointly with the second- and 
the fourth- boundary conditions in (|3.ip . we obtain that the modes v(-,-,k) are 
given by 

1 



v(t,x,0) =-&(*, 0) + (- 7 C/(^) 2 (t,0) - U(\nU)(Vy) 2 (t,0) -^(t,0))xe Ux 



v(t,x,0)=^<p t (t,0), x>0 
and 

v{t, x, k) =Ci tk e VkX + A k e Ux + B k xe Ux + C k xe UkX , x < 0, 

for k > 1, where 

(a + i)U 2 - U 

A k =i — (<Py) 2 {t,k) + aU$(t,k) + —(p t (t,k), 

Afe A k 

<yU 3 

Bk= 2 r ~ (Vy) 2 (t,k)+ 7 U 2 $(t,k), 



Ck = ^ Vk ^(t, k) + T f 2 " k 0(t, k), 

\ k {U-2v k )^ v> K l ' U-2uJ K 



° 2 - k ~\X k (U- 2v k ) + X k (U 2v k ) + K -U)(U- 2v k ) 2 ) Wy) [h 

jU 2 , jU 2 is k \ U{U-v k ) ^. 

V(*> k ) + 1 — 777 TT^^tit, k), 



U-2v k (C/-2^)V \ k (U-2v k ) 

Cl,k =C2,k — A k . 

The equation for the front now comes writing the last but one boundary condition 
in (13. ip . which we have no used so far, in Fourier variable, and taking advance of 
the formulas for the modes of u and v. It turns out that the equation for the front 
(in Fourier coordinates) reads: 

fo(«,0) + itf(S)a(t,0) = 0, 

U(U - v k ) ^ ( ~fU 2 !U 2 v k \ _ 

: ip f {t, k) + [ — — — + — — —j + v k -U\ (p(t, k) 



\ k {U-2v k y LX ' ' \U-2v k {U-2u k ) 

'^(JJ-Vk) jU 3 7 C/ 3 U(u k -U) 1 



A fe (C/-2^) \ k {U-2v k ) (v k - U)(U - 2v k ) 2 \ k 2 

or, even, in the much more compact form 

(X k U)fr(t,k) =^(U 2 -X 2 )(X 2 -jU 2 )ip(t,k) (3.3) 

+ \{Xl - WX 2 - ^U 2 X k + 4 7 C/ 3 )(^(t, k) 
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= (-4A2 + ( 7 -l)[/ 2 A fc )£(i,fc) 

+ ~ {Xl - 3UX 2 - 4 7 U 2 X k + 4 7 C/ 3 ) J^yf{t , k), (3.4) 

for any k > 0, if we set 

X k = ^U 2 + 4X k , k>0. 

Therefore, we have proved the following 

Proposition 3.1. Let (v,,v,<p) be a sufficiently smooth solution of System (|3.1[) 
such that the functions (x,y) i— > e~ Ux / 2 u{t,x,y) and (x,y) i— > e~ Ux ^ 2 v(t,x,y) are 
bounded in (— oo, 0] x [—£/2,£/2] and mix [—1/2,1/2], respectively. Then, the 
interface tp solves the equations (|3.4p /or any fc > 0. 

3.2. A fourth-order pseudo-differential equation for the front. Let us define 
the pseudo-differential operators (or Fourier multipliers) ££ and through their 
symbols, respectively 

b k = X k U, 

h = -4A£ + ( 7 - l)A fc t/ 2 , 
fk = \{X 3 k - 3UX 2 - 4 7 £/ 2 A fe + 4 7 £/ 3 ), 

for any k > 0. It is easy to see that 
SS = U(U 2 I ~4D yy )i, 

& = \{U 2 I - 4D yy )i - - A U{U 2 I - AD yy ) - jU 2 ^U 2 I-4D yy - U 
while the realization of J§? in L 2 is the operator 

L = -ADyyyy ~ ft - l^Dyy, 

with Hf as a domain. 

It follows from Proposition 13.11 that the front <p solves the equation 

±0(<p)=J?(<p) + &((<p y ) 2 ). (3.5) 

The main feature of Equation (13.51) is that the nonlinear part is rather unusual. 
Actually, it has a fourth-order leading term, as has. Therefore, (|3.5p is a fully 
nonlinear equation. More precisely, we have: 

Lemma 3.2. The operators SS and & admit bounded realization B : — > L 2 and 
F : H 3 — > L 2 , respectively. Moreover, B is invertible. 

Proof. A straightforward asymptotic analysis reveals that 
b k ~ 2^\~ k U, f k ~ 2X 3 k /2 , 

as k — > +oo, from which we deduce that and & admit bounded realizations 
B : Hj -> L 2 and F : Hf -> L 2 . 

Finally, since b k ^ for any fc > 0, it follows that B is invertible. □ 
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3.3. The third-order pseudo-differential equation for the front. In view of 
Lemma 13.21 we may rewrite Equation (|3.5p as: 

or, equivalently, as 

y t = sa?(y)+JZ{{v y ) 2 ). (3.6) 

We emphasize that Equation (|3.6[) is a pseudo-differential, fully nonlinear equation 
of the third-order, since the pseudo-differential operators stf and ^# have, respec- 
tively, symbol 

(U 2 - X 2 ) (X 2 - 7 [/ 2 ) Xl - 3UX 2 - ^U 2 X k + 4 7 C/ 3 

a k = 7777? j m k ~ 



WX k 1 AUX k 

for k > 0. Clearly, any smooth enough solution to (|3.5[) solves (I3.6[) as well. 
The following result is crucial for the rest of the paper. 

Theorem 3.3. The following properties are satisfied. 
(i) The realization A : L 2 of srf is a sectorial operator and the sequence 

(ttfe) constitutes its spectrum cr(A). In particular, is a simple eigenvalue of 
A and the spectral projection II associated with is given by 



Finally, a{A) \ {0} is contained in the left half-plane {A G C : RcA < 0} if 
and only if 7 < j c . 
(ii) The realization M : H 2 — > L? of the operator j% is bounded. 

Proof, (i). Let us split 

2A| -4A 2 fe + 4A 2 feV /£ + l-A fe [/ 2 + A fc7 [/ 2 2A J 

ak ~-— + uVTP+ir k : TT + 

for any k > 0. Since 

as fc — > +00, if 7 7^ 1/2, we can infer that the realization A of operator stf in i/j 3 

is well defined. Moreover, since A splits into the sum of two operators Aq (whose 

3 

symbol is (-2A| U^ 1 )) and A u which is a nice perturbation of Aq (being a bounded 
operator in Hi, which is an intermediate space of class Ji/ 3 between L 2 and H^), 
in view of [2UJ Prop. 2.4.1(f)] it is enough to prove that Aq is a sectorial operator. 
But this follows immediately from general abstract results (see e.g., [TTJ Chpt. 3]), 
or a direct computation. Indeed, if A has positive real part, then the equation 
Ait — Aqu = f has, for any / € L 2 , the unique solution 

+00 



u = R.(\,A )f = U^ }[k) 



t" XU + 2Xl /2 



and 



Proposition 2.1.1 in [3D] yields the sectoriality of Aq. 

Next we compute the spectrum of the operator A. Since is compactly embed- 
ded into L 2 , a (A) consists of eigenvalues only. We claim that a (A) consists of the 
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elements of the sequence (a*,). Indeed writing the eigenvalue equation in Fourier 
variable, we get the infinitely many equations 

A^(fc) - OfcVKJfc) =0, fc>0, (3.7) 

which should be satisfied by the pair A (the eigenvalue) and t/j (the eigenfunction) . 
It is clear that this system of infinitely many equations admits a non identically 
solution (ip(k)) vanishing if and only if A equals one of the elements of the sequence. 
The set equality a (A) — {a k ■ k > 0} is thus proved. 

Since the sequence (ofc) diverges to — oo as k — > +oo, all the eigenvalues of A 
are isolated. In particular, is isolated and, again from formula (|3.7[) . we easily see 
that the eigenspace associated with the eigenvalue A = is one-dimensional. To 
conclude that A = is simple, in view of [20j Props. A. 1.2 & A. 2.1], it suffices to 
prove that it is a simple pole of the resolvent operator. In such a case the associated 
spectral projection is the residual at A = of R(-, A). 

Clearly, for any A ^ o~(A), 

+oo 

for any (gl 2 . Hence, 

Ai?(A, A)( = 4>(0)w + V C(k)w k =: U( + Ri (A)C 

Since A ^ a k for any k > 1, and ak — > — oo as k — > +oo, there exists a neighborhood 
of A = in which the ratio |A|/|A — ak\ is bounded, uniformly with respect to 
k > 1. As a byproduct, in such a neighborhood of A = 0, the mapping A i-> i?i(A) 
is bounded with values in L(L 2 ). This shows that A = is a simple pole of operator 
A. 

To conclude the proof of point (i), let us determine the values of 7 such that 
a(A) \ {0} does not contain nonnegative elements. For this purpose, it suffices to 
observe that a k < for any k > 1 if and only if AX k + U 2 — jU 2 > for such k's, 
which is equivalent to 4Ai + U 2 — jU 2 > since (Afc) is nondecreasing sequence. 
Hence, the condition for a (A) \ {0} be contained in (—00, 0) is 7 < j c , where 

7c = 1 + ^- (3.8) 

(ii). As in the proof of Lemma (|3.2I) . it suffices to observe that mt ~ AfeC/ -1 as 
k -> +00. □ 

The linearized stability principle (see e.g., [3D1 Sect. 9.1.1]) and the results in 
Theorem 13.31 yield to the following stability analysis. 



Corollary 3.4. Let j c be given by (13.8[) . 

(a) If 7 < 7c, then the null solution to Equation (|3.6p is (orbitally) stable, with 
asymptotic phase, with respect to sufficiently smooth and small perturbations. 

(b) 1/7 > 7c, then the null solution to Equation (13.61) is unstable. 

4. Rigorous asymptotic derivation of the K-S equation 



The last question that we address is the link between (I3.6[) and Equation (|K-SI) . 
As in Section [2 we consider the small perturbation parameter e > defined by 
(see (|2TTjn 

7 = 1 + £. 
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Moreover, we perform the same change of dependent and independent variables as 
in (|2.2p . namely: 

, T r) e 

The key-idea is to link the small positive parameter e and the width of the strip, 
which will blow up as e — > 0. For £ > fixed, we take I of the form: 

£ " y/eU' 

hence 7 C (see (|3.8[1 ) converges to 1 as e — > 0. 

In view of Corollary 13.41 in order to avoid a trivial dynamics, we assume that 
7 C > 1. This means that we take the bifurcation parameter £q larger than Air and 
obtain that j c G (1, 1 + e). 

In the new variables is replaced by the operator S$ z = U 2 ^J(I — AeD, m ). 
Lemma 13.21 applies to this operator and guarantees that, for any fixed e > 0, the 
realization B e : — > L 2 of SB e is bounded. However, the perturbation is clearly 
singular as e — > 0, since obviously B e — > U 2 I. Therefore, it is hopeless to take the 
limit e — > in the third-order equation (|3.6p . whose behaviour is clearly singular 
as e — > 0. Fortunately, the fourth-order equation (|3.5[) is more friendly, since, after 
division by e 3 and U 3 , it comes: 
Q f x 

y/1 - AeD m V = -iDrjrjrjrjIp - D^t/j 



Ot 

+ &D m )i - 3(1 - 4sD vv ) - 4(1 + e) [y/I - AeD m - ij } (D^) 2 , 

(4.1) 

which is the perturbed equation we are going to study, with periodic boundary 
conditions at r) = ±£q/2. 

Mimicking (|3.5I) . we rewrite (|4.ip in the abstract way: 

^-^ = ^ + J^((^) 2 ), (4.2) 
where the symbols of the operators 3§ e , Jz? and are respectively: 

b s ,k = Xe,ki 

s k = -Afc(4A fe - 1), 
1 



f e , k = -{Xl k - ZX 2 Sik - 4(1 + e)X e , k + 4 + 4e), 



for any k > and 



X e , fe = v / l + 4eA fc , k > 
Writing (|4.2p in discrete Fourier variable gives the infinitely many equations 

b £ ,k$r{T, k) = -Afe(4A fe - l)^(r, k) + f e ,kl^(T, k), 

for any k > 0. Note that the leading terms (namely at order in e) of b e ^k and / Ei fe 
are 1 and —1/2, respectively. 

Fix T > 0. For $ € HJ™ (m > 4) the Cauchy problem 

*t(t,j?) =-4$ ww (r,?7)-$ w (r,r?)- i($ I( (r, t?)) 2 , t > 0, |r?| < ^, 

< D*#(r,-V2)=^(r,V2) ) r > 0, fc > 0, 

t $(0,7?) = $0(77), |r?|<f. 



ASYMPTOTIC ANALYSIS IN A GAS-SOLID COMBUSTION MODEL 



13 



admits a unique solution $ <E C([0, T]; Hf) such that $ T G C{[0,T}; H™~ 4 ) (see 
e.g., App. B]). 

Through <J> we split ip = $ + ep e - For simplicity, we take zero as the initial 
condition for p e and, to avoid cumbersome notation, in the sequel we usely write p 
for p e . 

If i/} solves (pl~2]) . then 

J^e(p) + -Wr) = JSf(p) + -C(($„) 2 ) + £^e((p I( ) 2 ) + ^e(%P V ), (4.3) 

where the symbols of the operators J4? £ and M e are 

k.fc = \( x z,k - 1), m E>fe = ^(Xl k - 3A £ 2 fc - 4(1 + e)X £tk + 6 + 4e), 
for any k > 0. 

Proposition 4.1. There exists a positive constant C* suc/i £/ia£ £/ie following prop- 
erties are satisfied for any e £ (0, 1]: 

(a) for any s = 2,3, . . ., the operators SS e and Jif e admit bounded realizations B e 
and H E , respectively, mapping into H^~ 2 . Moreover, 

W B z\\l(h$,h%- 2 ) + W h cWl(h> v h;- 2 ) ^ c *- 

Finally, the operator B e is invertible from HS to H^~ 2 ; 

(b) for any s = 3,4,..., the operators J? e and ^ e admit bounded realizations F e 
and M e , respectively, mapping H s intoH s ~ 3 . Moreover, 

\\ F e\\L(Hf,Hf- s ) + ll M elli(fl-jf,Hj- 3 ) - C *' 

Proof. The statement follows from an analysis of the symbols of the operators SS e , 
j£" e , ffl e and jM z . Without much effort one can show that 

|/i e ,fc| < 4A fc , |m E , fe | < 2A| + 25A fc , 

for any k > and any e G (0,1]. These estimates combined with the formulas 
7^ b E: k = eh e ,k + 1 and fk = em e ^ — 1/2, for any k > and any e G (0, 1], yield 
the assertion. □ 



Instead of studying Equation (13.51) , we find it much more convenient to deal with 
the equation satisfied by £ := p„, i.e., 

J^e(C) + *%{*t) = J2?(C) + ^((* 2 ) I; ) + £^((C 2 ),) + 2^ e ((*C)v). ( 4 - 4 ) 
which we couple with the initial condition £(0, •) = 0. Here, W = 

4.1. A priori estimates. For any n = 0, 1, 2, . . . and any T > 0, we set 

X„(T) = {C G C([0,T];ffp)nC 1 ([0,T];I 2 ) : ( T € C([0, T]; H„ 2v(n+1) )} , 



For any e > 0, we introduce in hV 2 the norm ||£||f = J2k=o + 4eAfc |C(fc)| 5 



where a V 6 := max{a, 6} 

Arln/in in / / ^ /*" +- Via n /-A v m 1 1 /" 1 1 2 \ 

;k=0 
1/2 

for any £ G HJ . Note that, for any fixed e > 0, || • ||i s is a norm, equivalent to 

1 /2 

the usual norm in H J . 

The main result of this subsection is contained in the following theorem, where 
we set * = ($oV 
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Theorem 4.2. Fix an integer n > and T > 0. Further, suppose that 'J'o € . 
Then, there exist E\ = E\{n,T) G (0,1) and = K n (n,T) > suc/j f/ia£, £/ 
C € A„(Ti) is a solution on the time interval [0,Ti] of Equation (|4.4|) /or some 
Ti < T, £/ien 



whenever < s < e±. 

Note that the assumptions on guarantee that 'J belongs to C([0,T]; H? +4 ) n 



The proof of Theorem 14.21 heavily relies on the following lemma. 

Lemma 4.3. Let Aq, cq, c\, ci, C3, e, Tq, Ti be positive constants with T\ < Tq. 
Further, let f £ ,A £ : [0,Ti] — > R be a positive continuous function and a positive 
continuously differ entiable function such that 



Then, there exist E\ = si(Tq) G (0, 1) and a constant K — K(Tq) such that A £ (t) < 
K for any t € [0,Ti] and any e € (0,£i]. 

Proof. When f £ identically vanishes, the proof follows from [2, Lemma 3.1], which 
shows that we can take £i(T ) = 3c^/(16cic 3 (e C2T ° - 1)) and K < 4cie C2T °/(3c 2 ). 

Let us now consider the general case when f e does not identically vanish in [0, T{\. 
We fix £ = £q(Tq) < 3c^/(16cic 3 (e C2T ° - 1)) such that 9c c^ - l2 Cl c 2 e C2To e - 
16cfe 2c2T °e§ > 0, and e G (0,e ]. We claim that c - s 2 (A £ (t)) 2 > for any 
i-e[0,Ti]. 

Let (0, T e ) be the largest interval (possibly depending on e) where Co — sA £ — 
£ 2 {A £ ) 2 is positive. The existence of this interval is clear since A £ vanishes at 0. 
The positivity of c — e 2 (A £ ) 2 in (0, T £ ) shows that A' £ < c± + c 2 A £ + c 3 sA 2 in such 
an interval and, from the above result, we can infer that A £ (t) < (4cie C2T °)/(3c2) 
for any r £ [0, T £ ], so that c - s 2 (A £ (T £ )) 2 > 0. By the definition of T £ , this clearly 
implies that T e = T x . □ 

Proof of Theorem 14.21 Throughout the proof we assume that T± < T is fixed and 
e and t are arbitrarily fixed in (0, 1] and in [0, Ti], respectively. Moreover, to avoid 
cumbersome notation, we denote by c almost all the constants appearing in the 
estimates. Hence, the exact value of c may change from line to line, but we do not 
need to follow the constants throughout the estimates. We just need to stress how 
the estimates depend on e. As a matter of fact, all the c's are independent not only 
of e but also of r, ^> and C- On the contrary, they may depend on n (and, actually, 
in most the cases they do). Finally, we denote by K^) a constant which may 
depend on n and also on As above K(fy) may vary from estimate to estimate. 

The first step of the proof consists in multiplying both sides of Equation (|4.4[) by 
(— l) n D 2n ( and integrating by parts over (~l /2, to/2). This yields to the equation 



sup ||D^C(r,-)lli B <^i 



(4.5) 



G 1 



([0,T];7J"+ 2 ). 





2 



(4.6) 




(H e (? r (r, •)) - M £ ((* 2 Ut, •))) (-1)"-Df C(r, -)dn 
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in 



+ e I F £ ((Ov(T,-))(-l) n D?aT,-)dr) 



+ 2 F E ((*C),(T,-))(-l) n J DfC(r,-)^ 



(4.7) 



(4.8) 



Using Parseval's formula and the definition of the symbol b St k> one can easily 
show that 



1 d 



B £ (C T (r, ■))(-l)"< l C(r, -)dn = - — II^CCt, OH* 



(4.9) 



We now deal with the other terms in (|4.8|) . Integrating n-times by parts and, 
then, using Poincare-Wirtinger and Cauchy-Schwarz inequalities, jointly with Propo- 
sition HUl it is not difficult to show that 



<*r(*) + |£>«c(T,oii, 

(4-10) 



(fle(* T (r, r,)) - M e ((* 2 )„(r, •))) (-l) n ^«C(r, 0*7 



for any ( € X n (T x ). 

Estimating the other two integral terms in the right-hand side of (|4.8[> demands 
some more effort. The starting point is the following estimate 



F £ ( Xv (T,-))(-lTD 2 n n C(T,-)dv 



<cei\D:; + \(r, -)|2K +2 C(t, 0|2 + cs\D-+ 2 x(r, -)| 2 |^ +1 C(^ Ola 

+ cV-e\D- X (r, 0| 2 |^ +2 C(r, Ola + c\D^ X (r, OH^Cfr Ob, (4.11) 

which holds true for any X G C([0, 7\]; i7j 4v2n ). Such a formula follows observing 
that 

/ o 

n T-\2n / 



F e { Xv {T,-)){-l) n DfQ{T,-)dr, 



-HOC: 



<£ASi/ e ,*ii^M)iic(r,AOi 

k=0 

<c£ \ n k {si XI + e\ k + e*A* + l) |^(r, fc)||C(r, fc)|, 

fc=0 

and, then, using Young inequality, to estimate the terms in the round brackets. 
Now, we plug x — ( 2 m t° (|4. and use the estimates 

\D^ +2 (C(r, ■))% < C (lA" +2 C(r, Ola|^C(r, Ola + lA" +1 C(r, Ola), (4-12) 

\D^C{t,-)Y\ 2 <c\D-C{t,-)\1 (4.13) 

(which can be obtained using the Poincare-Wirtinger inequality and Leibniz for- 
mula), the Cauchy-Schwarz inequality and, again, the Poincare-Wirtinger inequal- 
ity, to obtain 



F £ ((C 2 WT,0)(-l)^ n C(T,0dr? 



<ce-*\D"((T, OlaPJTCOr, Oil + ce^D^^r, OlaW^r, Ob 
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+ ce\D-C v (r, ■)\2\D^\(t, 0| 2 K +2 C(t, Ola + ^ +1 C(r, .)\% 

+ cei (1 + e)\D»ar, OIII^CCt, Ola + c\D^(r, OIll^C^ Ola 
<cei|^C(r, -)|2|^ +2 C(r, Oil + <** \D-+\(t, -)|H^ +2 C(^ Ola 

+ c £ |^ +1 C(r, OI2I A7 +2 C(r, 01= + ^ +1 C(t, 01 WC(t, Ola 

+ ce* |^C(r, 0III^ +2 C(r, Ob + c|D»C(r, Olll^Cfr -)b 
<cet |^C(r, 0| 2 |^ +2 C(r, Oil + ce 2 \D^\(r, 0| 2 + cs\D^((t, -)\% 

+ C |^C(r,0|^ + C |^; +2 C(r,OI^ 

<cei \D;C(r, 0| 2 K +2 C(t, Oil + ce 2 \D-+\(r, 0| 2 + c\D^+ 2 ((t, Oil + c|^C(r, 

(4.14) 

In the similar way, using the estimate |-D™(#C) 0, |a < c|£>™C(t, |a \D™V(t, ■) |a, 
(with m £ {n, n + 2}) in place of (I4.12|) and (14 131) . from (|4.11[) we get 

J ^F £ (^CUrr))(-l) n D^C(r,-)d V 

<cei \D n n +H(r, 0| 2 |^ +2 C(r, Oil + cs\D^+H(r, 0|a|^ +1 C(r, Oil 

+ c £ |^ +2 *(r, 0| 2 |^ +2 C(r, OH + c\D^{t, 0| 2 |^C(r, Oil 

+ celX^*^, 0| 3 |£>? +2 C(t, Oil + cff- 1 !^*^, OIII-O^CCt-, Oil + cff|Z^+ a C(^r, Ola. 

(4.15) 

for any 6 > 0. We just mention the inequality 

I^^OIap^Olap^C^OIa 
KcS-^^t, -)\ 2 2 \D n C(r, 0| 2 + S\D n+1 C(r, Oil, 

obtained by means of Young and Poincare-Wirtinger inequalities, which we use to 
estimate one of the intermediate terms appearing in the proof of (|4.15l) . 
Now, taking c6 — 5/2, we get the estimate 

<K(V) (e\D^((T,-)\l +e\Di;+ 1 ar,-)\l + |^C(r,0ll) + ||^ +2 C(^,0il- 

(4.16) 

From dUSD , (|4~§1) , (ETTUj) . P~14"l) . (|4~I5|) and the interpolative inequality 

l^ +1 C(r,0ll < |d»c(t,0|1 + ^™+ 2 c(t,0|1, 

we can infer that 

~KC(t, 0II|, £ + (1 - (*) - cellar, 0|a)|^ +2 C(r, Oil 

+ A"(*)|D»C(r, Oil + £ if(*)|^ +1 C(r, Oil 
+ C£ |^C(r,0l^ + c £ 3 |^ +1 C(r,0l^ 
<K(9) + K(*)\DWt, OH + (*)|^ +2 C(r, Oil 
+ ce|^C(r, Oil + c£ 3 |^C(r, Olal^ +2 C(r, Oil, 



ASYMPTOTIC ANALYSIS IN A GAS-SOLID COMBUSTION MODEL 



17 



which we can rewrite in the form 

^IKC(r, OII| lB + (2 - eK(V) ce 2 \\D-C(r, .)||| i(S ) \D^ 2 ((r, 

+ K{*)\\mC(r, + cellar, -)ll i, e . 

estimating 2\\D%((t, -)\\i jE < \\D^(t, -)ll i £ + 1 and recalling that e € (0, 1]. 

Up to replacing (0, 1] by a smaller interval (0, £0], we can assume that eif(^) < 1 
for any e € (0, £0]. Hence, applying Lemma 14.31 with 

Co = l, Cl C 2 c 3 = c, 

A e (r) = ||^C(r,-)ll|, e , / E (r) = K +2 C(r,-)li 
we conclude the proof. □ 

Now, taking advantage of the previous a priori estimates, which can be extended 
also to variational solutions to (I4.4p belonging to the space spanned by the 
functions w\ , . . . , wn (with constants independent of N G N), and using the classical 
Faedo-Galerkin method, the following result can be proved. 

Theorem 4.4. Fix T > 0. Then, there exists £o{T) > such that, for any < 
£ < £o(T), Equation (|4.4p /ias a unique classical solution £ on [0, T], vanishing at 

T = 0. 

4.2. Proof of Main Theorem. Since the unique solution £ of Problem (|4.4j) is 
the candidate to be the 77-derivative of the solution p to Problem (|4.3[) . p should 
split into the sum p(r, rj) = (3^(C,))(t, rf] + v(t) for some scalar valued function v, 
where 

1 / 2, 



(^(C))(r, r?) = y ^ C(s)ds - - C(«) - Tq ) da. 

Imposing that p in the previous form is a solution to (|4.3p and projecting along 
LT(L 2 ), we see that p is a solution to (|4.3[) if and only if u solves the following 
Cauchy problem: 

^ = _n(if £ (^))-i £ n(c 2 )-n($, 7 c) ) 

u(0) = 0. 

Since this problem has in fact a unique solution, and ^(C) + u vanishes at r = 0, 
we conclude that problem (|4.3p is uniquely solvable. 

To complete the proof, we should show that there exists M > such that 

sup \ P {t, V )\<M, (4.17) 

rE[0,T] 

v e[-e /2,e /2] 

uniformly in < e < eo(T). Once this estimate is proved, coming back from 
Problem (|4.3[) to Equation (jl.lip . we see that the latter one has a unique classical 
solution ip : [0, JF^ } x K — > R, which is periodic (with respect to the spatial variable) 
with period £ e = £ /(^/eU), and satisfies </j(0, •) = eJ7~ 1 $ (v / ^C^')i as weu as the 
estimate 

||^(i,-)-et/- 1 $(te 2 C/ 2 ,-ViC/)||c(M e /2A/2]) < — te[0,T e ], 
as it is claimed. 

So, let us prove (|4.17p . For this purpose it is enough to use the a priori estimate 
(|4.5p jointly with the Poincare-Wirtinger inequality, to estimate f, and just (I4.5[) 
to estimate u. This completes the proof of the Main Theorem. 
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5. Numerical experiments 

In this section, we intend to solve numerically Equation (I4.1j) for small positive 
e and illustrate the convergence to the solution of K-S equation. 

In order to reformulate (I4.1j) on the interval [0, 2tt) with periodic boundary con- 
ditions, we set x — r]/(2£o), where £ — Iq/Att. It comes: 



d e \ 1 1 



Ot 



to / 4£ U 



I-4^D XX \ -3(7- 4iDj\ -4(1 + e){ I--^D xx -l) 



io J V *o J VV £ o /J 164 

Next, we define the bifurcation parameter j3 = 4£o as in [13], [7]. After multipli- 
cation by /3 2 , it comes: 

d 



— [^^ - 4ep3D xx ) i> = -AD XXXX ^ - f3D xx i; 



1 1(7 - ^7> xa )S - 3 (7- ^23 as J - 4(1 + e) ^1 jD xx -I}} ID, 



Finally, we rescale the time, setting t — t/(3 



2. 



d ( / 4e . 

g- t ^-7> KX ) V = -4D xxxx ip - PD xx ip 



+ \ I (7 - ^A«J ' - 3 (l - |7J xa J - 4(1 + e ) (Jl-^ Dxx -I)} (D X ^Y 

and setting e' = e//3, the prime being omitted hereafter, we obtain: 
d 

> 



— yyi - AeD xx j if) = -4:D xxxx ip - /3D xx ip 



+ J {(7 - 4e7> x:c )i - 3(7 - 4eD xx ) - 4(1 + e) (yi-4eD xx - 7) } (7>^) 2 . 

(5.1) 

The initial condition is given by ij)(0, •) = ipO: where i/'o is periodic with period 2tt. 
Note that, in contrast to [T3], [7], we do not subtract the drift. 
Equation (|5.1I) reads in discrete Fourier variable: 

d 



— (V 1 + 4efc 2 ) #(t, fc) = -4fc 4 ^(i, k) + (3k 2 $(t, k) 



+ I {(1 + 4efc 2 )t - 3(1 + 4 £ A: 2 ) - 4(1 + e) (V 1 + 4efc 2 - l) } (V>*) 2 (t, k). 

We use a backward-Euler schema for the first-order time derivative, treat implicitly 
all the linear terms and explicitly the nonlinear terms. The implicit treatment 
of the fourth- and second-order terms reduces the stability constraint, while the 
explicit treatment of the nonlinear terms avoids the expensive process of solving 
nonlinear equations at each time step. For simplicity, in the rest of this section we 
use the notation f k instead of f(k). It comes: 

(7l + 4efc 2 ) rk+ ' At ^ = -4fc% +1 + (3k 2 r k +1 



+ f {(1 + 4efc 2 )^ - 3(1 + 4efc 2 ) - 4(1 + E ) (V 1 + 4 £ £: 2 - l) } {[(^)'f } k , 
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where {(ip x ) 2 }k represents the fc-th Fourier coefficient of (4>x) 2 - This method is 
of the first order with respect to time. From the previous equation it is easy to 
compute the A;-th Fourier coefficient ■ One gets: 

x ((1 +4efc 2 )^™ 

+ ^|(l + 4eA: 2 )i - 3(1 + 4e/fc 2 ) - 4(1 + e)[(l + 4efc 2 )^ -1]} 

x{[(^) n ] 2 } fc ), (5.2) 

Practical calculations hold in the spectral space. We use an additional FFT to 
recover the physical nodal values ipj from tj)^, where j stands for the division node 
in the physical space. 

The numerical tests aim at checking the behavior of the solutions of Equation 
(|5.1[) for values of e close to 0, and compare them to the Kuramoto-Sivashinsky 
equation. In Figures OH and [SJSJ we plot consecutive front positions computed 
using Equation f|5.2[) taking /? = 10, 20 and giving to e the following values: 0.1, 
0.01, 0.001, and (which corresponds to the K-S equation). 




FIGURE 1. Front propagation with j3 — 10, e — 0.1 and ipo(x) = sin(x). 




FIGURE 2. Front propagation with (3 = 10, e = 0.01 and ipo(x) = sin(a;). 




Figure 3. Front propagation with f3 = 1Q, e = 0.001 and tpo(x) = sin(a^). 
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Figure 4. Front propagation with (5 = 10, e = and Vo^) = sin(x). 
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Figure 5. Front propagation with /3 — 20, e = 0.1 and i/>o(x) = sin(x). 
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FIGURE 6. Front propagation with (3 = 20, £ = 0.01 and ^o(^) = sin(x). 




Figure 7. Front propagation with = 20, £ = 0.001 and Vo(#) = sin(x). 
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Figure 8. Front propagation with (3 — 20, e — and ipo(x) — sin(x). 

We now investigate the dynamics of Equation (I5.1[) with respect to the parameter 
(3. For this purpose, we fix e = 0.001. 

The numerical simulations confirm that, as for the K-S equation, turns out 
to be a global attractor for the solution to Equation (|5.ip . for any (3 £ [1,4]. A 
non-trivial attractor is expected for larger /3's. In Figures |9"1TT21 we can see the front 
evolutions generated by (I5.2j) with (3 = 30, 60 for two different initial conditions. 
In all the below figures, the periodic orbit is clearly observed. 




Figure 9. Front propagation with /3 = 30, e — 0.001 and ipo(x) = sin(a;). 
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FIGURE 10. Front propagation with j3 = 30, e — 0.001 and 
ipo(x) = cos(x). 
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Figure 11. Front propagation with j3 = 60, e = 0.001 and 
ipo{ x ) = sin(x). 
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Figure 12. Front propagation with j5 = 60, £ = 0.001, and 
ipo(x) = cos(x). 

Summing up, our numerical tests confirm that Equation (15. ip preserves the same 
structure as K-S equation. Larger (3 generates an even richer dynamics, see Figure 
1131 where the front propagation is captured from a computation with /3 = 108. 
As predicted in [13], the front evolves toward an essentially quadrimodal global 
attractor. 




Figure 13. Front propagation with j3 = 108, e = 0.0001, and 
tpo(x) = 0.1(cos(cc) + cos(2ir) + cos(3a;)). 
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